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Abstract. In this work we analyze the relation between the multiplicative 
decomposition F = F e F p of the deformation gradient as a product of the 
elastic and plastic factors (12 17 , 18 , 26 ) and the theory of uniform materials 
( |33l I40l I41 p , We prove that postulating such a decomposition is equivalent 
to having a uniform material model with two configurations - total <j> and the 
inelastic <pi. 

We introduce strain tensors characterizing different types of evolutions of 
the material and discuss the form of the internal energy and that of the dissi- 
pative potential. The evolution equations are obtained for the configurations 
((f>,<j>l) and the material metric g. 

Finally the dissipative inequality for the materials of this type is presented. 
It is shown that the conditions of positivity of the internal dissipation terms 
related to the processes of plastic and metric evolution provide the anisotropic 
yield criteria. 
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1. Introduction. 

The objective of this work is to investigate the relation between the geometrical 
theory of uniform materials and the multiplicative elasto-plastic decomposition of 
the deformation gradient of Bilby-Kroner-Lee (BKL-decomposition) and Nemat- 
Nasser (see (2j [HI [HI [26]). 

Such a relation was first studied in [29] . In particular, the relation between 
the inhomogeneity velocity gradient L p (see below) and the plastic distortion rate 
L = F p • (F p ) _1 was introduced. In this paper we study the geometrical form of 
the relation introduced in |29j . 

In Section 2 we introduce the basic concepts and review properties of uniform 
materials. In Section 3 a bijective correspondence between the BKL decompositions 
of the gradient of a configuration </> of an elasto-plastic solid and the triple (<f>, cf>i, P) 
is established. Here P represents the uniform material structure and 4> and </>i 
are respectively, total and inelastic (intermediate) material configurations. 

In Section 4 we introduce the natural strain tensors measuring the relations 
between the Cauchy-Green deformation tensors C(4>) and C(cf>i) and the material 
metric g induced by the uniform structure P. In the same section the combinations 
of these tensors, material metric and its curvature characteristic independent on the 
decomposition of plastic deformation gradient F p — o D are determined and 
the strain rate tensors are introduced. 

In Section 5 the form of internal energy u depending on variables (4>, cf>i, g) and 
their derivatives is postulated and the dissipative potential V is introduced. We 
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also formulate the system of equations describing evolution of dynamical variables 
((f), (f>i,g,). In the same section different stress tensors present in our scheme are 
defined and relations between them arc discussed. 

In section 6 we write down the dissipative inequality for the suggested scheme 
and separate the terms corresponding to the internal dissipation related to the 
processes of integrable inelastic and uniform structure evolutions. We show that 
the conditions of positivity of the corresponding terms in dissipative inequality 
provide the anisotropic yield criteria for initiating the corresponding processes. 

Another form of a relation between the finite elasto-plasticity based on the mul- 
tiplicative decomposition and the uniformity structures using the second-order con- 
nection was suggested by S. Cleja-Tigoui, see [7]. 

2. Uniform Materials: material connections and material metrics. 

Uniform materials enter the scene of material science about 1952 when K.Kondo 
introduced the material connection and the material metric as the tools to model a 
properties of materials. Later development in the works by K.Kondo, B.Bilby and 
his collaborators, W.Kroner, W.Noll ([33]) and C.C. Wang (see @0j|4T]) establish 
the basis of this theory. In the works of 1980-present by M.Elzanowski,M. Epstein, 
M. De Leon, G. Maugin different aspects of this theory: models of higher grade 
uniform materials, dynamics of material properties, thermodynamical properties 
of such materials, role of Eshelby stress tensor, geometry of functionally graduate 
material, etc., were further developed. 

In this Section we present the basic geometrical structures of the theory of uni- 
form materials that will be used in latter parts of the paper. Our presentation is 
based on El QJ [Ml- 

2.1. Material and physical spaces. A material body (material manifold) is 

usually represented by a connected 3-dimensional smooth oriented manifold M with 
a piecewise smooth boundary dM . Constructions of this paper are local, so it is 
sufficient to consider M as a connected open domain in R 3 with local coordinates 
X 1 , 7=1,2,3. 

As the physical space our body is placed in we consider the 3-dimensional Eu- 
clidean vector space (E 3 ,h), h being the (flat) Euclidean metric. We introduce a 
global Cartesian coordinates x l in M 3 . In these coordinates the metric h takes the 
form h = hijdx l dx 3 . 

We will also use the concept of "archetype" ([9], [IS]), a 3-dimensional vector 
space V endowed with a standard Euclidean metric and the orthonormal basis 
eo = {ei, i = 1,2,3}. For convenience we identify the "archetype" space V (see 
[12"ll29j) with the tangent space at the origin O of the physical space: V = To(M 3 ) 
and its metric with the metric h at the origin. 

2.2. Configurations and the Cauchy metric. Configuration of the body M 
is a (diffeomorphic) embedding <fr : M — > E 3 into the physical space E 3 . To each 
configuration <f> there corresponds the deformation gradient - the mapping from 
the tangent space Ty(M) at the point X 6 M to the tangent space T^ X )(E 3 ) at 
the point 4>(X) S E 3 , [22], 

F(X) = fax : T X (M) T HX) (E 3 ), 
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given, in coordinates X A , x % , by the matrix of partial derivatives 

F(X)} = 4> l j- 

Here and below we will use notation tffj = ■^y(X) for the partial derivatives of 
configuration components (j) l {X). 

To a configuration 4>{X) there corresponds the right Cauchy-Green deformation 
tensor - the flat metric C(4>) = (f)*h in M obtained as the pullback of Euclidian 
metric h in physical space by the configuration mapping cf>. In coordinates (X 1 ) 
tensor C((f>) has the form 

CWv = h.X' ,<:<[,. (2.1) 

We will fix a specific configuration </> Q and call it the reference configuration. 
Usually it presents the state of the material body that is free from loads and stresses 
(see [39l [25]), although it might happen that such a configuration does not exist 
and one has to choose a reference configuration differently. The body M is often 
identified with its image under the embedding <j> . 

To the reference configuration <j) there corresponds its Cauchy-Green tensor 
called the reference metric in M: 

go = C(<t> a ), g oIJ = hij^^j, (2.2) 

and the corresponding reference volume form v a (X) — y/~\g ~\dX 1 A ... A dX n . 

Using the mapping inverse to the reference configuration 4> : M — > E 3 one can 
define the frame p Q in M by the rule 

Po{X) = ^ x (e ), (Poh = qx* ' 1 = 1;2 ' 3 - 

Q 

From now on we assume that the coordinates X 1 are introduced in the material 
manifold M using the reference configuration, i.e. X I (X) — (f>l(X). Then the 
vectors of the frame p take the form (p )i — g^Ti ^=1)2, 3. 

Finally we define a history of deformation as a time parameterized family of 
smooth configurations: 4>{t, X) : M x E — > E 3 . 

2.3. Uniform materials, I. Recall ([40[[33]) that a material is called hyperelas- 
tic if its constitutive response (to a loading conditions) at any configuration <f> is 
completely characterized by two scalar functions: 

(1) The elastic energy density function (per unit of reference volume v ) 
W(X,F(X)) depending on a material point X e M and the deformation 
gradient F(X) at this point; and 

(2) The mass density function p re f(X) > in the reference configuration <f> . 

Next we introduce the basic notion of a uniform material (body). Intuitively 
speaking, a uniform body is one that is made of the same material at all its 
points. The property of uniformity is characterized in terms of a parallelism K\ in 
the body M ([40l [39l [9]). More specifically, a hyperelastic material body (M, W) is 



Here and bellow for a diffcrentiablc mapping tp : M —> N between manifolds M and N we 
denote by ifi*x '■ Tx(M) — » T^(x){N) the linear mapping of tangent spaces at a point X £ M. 
In coordinates [X 1 ,x l ) mapping ip*x is given by the matrix Fj = (jf*. Corresponding mapping 
of the tangent bundles will be denoted by V*: 4>* ■ T(M) -* T{N), see [22], Ch.l.,SMK 
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called uniform if for any two material points X, Y there exists a linear isomorphism 
K\ : Tx(M) — ► Ty(M) between tangent spaces at these points such that 

K x *(W(Y,F(Y))dv (Y)) = W(X,F(Y)oK%)dv (X) (2.3) 

for all values of deformation gradients F(Y) at Y. Here K x * is the pullback of the 
n-form of energy density by the mapping K x . 

Introduce the scalar factor X x , characterizing the behavior of the reference vol- 
ume form under the parallelism K x : K x *v (Y) — X x v (X). Then, in terms of the 
energy density function W condition (2.3) takes the form 

X x W(Y,F(Y)) = W(X,F(Y)oK x ) (2.4) 

for all points X, Y in M and for all values of deformation gradient F(Y) at the 
point Y. 

2.4. Material connections. The localization of the definition of uniform mate- 
rials given above leads to the introduction of a linear connection (material con- 
nection) to in M having vanishing curvature (an absolute parallelism, see ([16J)). 
Having such a connection available, the mappings K x are defined by the parallel 
translation defined by connection to from the point X to the point Y along any 
curve connecting X and Y (result of such translation is independent on the choice 
of a curve due to the vanishing of the curvature). The torsion tensor T of connection 
to provides the measure of non-homogeneity of the material, see (5J|5]. 

It is known (see [16], Ch.2) that in a simply connected body M which admits 
a global tangent frame, a zero curvature connection is determined by a choice of a 
global tangent frame parallel with respect to the connection to 

p(X) = { Pk = P [{X)d XI , k = 1, . . . , 3, V" Pfe = 0}. 

Remark 1. A choice of such a frame is unique up to the (natural) right action of 
the group GL(n, R) on the tangent frames and the left action of the symmetry 
gauge group G M of the connection to (see ([HUM]))- 

A global frame p may also be defined by the uniformity mapping smoothly 
depending on the point X 

Px-V^Tx(M), Px{e i = (Px) I i d x i, i = 1,2,3. (2.5) 

Mapping P x defines the linear isomorphism of the archetype space V with the 
tangent space at each point X € M . Section p and the uniformity map P are 
related by 

p(X)=P x (e )^ Pl (X)=Pf-^j. (2.6) 

Parallel translation K\ defined by the connection u> can be written in terms of 
the uniformity mapping as the composition 

K Y x =PyoP x \ 

Using the reference frame p (see above) and the frame {ej} in the space V, one 
can associate to a material frame p two other geometrical objects: 

(1) A smooth mapping k : M — > GL(V),X k(X) (an element of the gauge 
group GL(V) M ) such that for all X G M 

P./(A) = Po j(X) ■ k(X) & p*j(X) - (pf jk(X)i, I, J = 1, 2, 3, 
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here GL(V) is the group of invertible linear transformations of the archetype 
space V; 

(2) A non-degenerate (l,l)-tensor field Dj(X) such that 

D(A) Po (A) = p(X),i.e. p((X) = D(X)5(p )f PO = Df(X), i, I = 1, . . . , 3, 

last equality being true due to (po)l(X) = Sf . 
Non-degeneracy of the (l,l)-tensor D(X) means that D(A) G GL{T X {M)). 

Using the relation between the frame p and the corresponding gauge mapping 
k : M — > GL{V) we get the relation between k and the uniformity mapping P cor- 
responding to the frame p, namely, Pi(X) = Px{&i) = (po)ik(X) = P ,x(ei)k(X), 
so that 

Px =Po,xok(X). 
These considerations are summarized in the following 

Proposition 1. Let M be a simply connected parallelizable (i.e. admitting a global 
frame ) manifold . With a choice of a reference configuration (f> and a frame in 
the archetype space V there is a bijection between the following objects: 

(1) Global frames p in M (global smooth sections of the frame bundle F{M) ); 

(2) Smooth uniformity mappings Px ■ V — > Tx(M); 

(3) Smooth mappings k : M — > GL(V),X — » k(X) (elements of the gauge 
group GL(V) M ) such that for all X e M 

p(X) - p (X)k(X); 

(4) Non- degenerate smooth (1,1) -tensor fields Dj(X) in M such that 

T>(X) Po (X) = p(X), 
or, in terms of uniformity mappings P and P 
B(X) = P x oP- 1 . 

Remark 2. It is the bijection between the first two and the last types of geometrical 
objects (non-degenerate (l,l)-tensor fields) that will be primarily used in this paper. 

2.5. Uniform materials, II. A uniformity mapping P determines its own vol- 
ume form by translating to the material the Euclidian volume element from the 
archetype: vp(X) = P^ 1 *(e\ Ae2 A e3). Denote by Jp(X) the factor relating two 
volume forms v and vp 

v P (X) = J P (X)v (X)- 

-Jacobian of the mapping 

Comparing definition of the factor X x in (2.4) with the definition of the factor 
Jp(X) we get, for a uniform material following relation between these factors: 

Y _ Jp(X) 

Xx - JpIy)- (2J) 

In terms of the volume factor fip uniformity condition (12. 4|) takes the form 

J P 1 (Y)W(Y,F(Y)) = J- 1 (X)W(X,F(Y)oP Y o P^) (2.8) 

Combining the deformation gradient F(X) and the uniformity mapping Px one 
gets the linear automorphism of the archetype space Ax = F(X) o P x £ GL(V). 
Comparing (2.4) with (|2.8p we rewrite the condition (|2.8p as follows 

Jp\X)W(X,F(Y)oP Y oP x 1 ) = J p 1 (Y)W(Y,F(Y)) = J P 1 (Y)W(Y, F(Y)oP(Y)oP(Y)- 
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for arbitrary points X, Y £ M and an arbitrary value of the deformation gradient 
F(Y) at the point Y. 

Define a function W of a point X S M and a linear mapping A G G-L(V) by 
setting 

#(JT,A) = J P 1 (X)iy(X,AoP x 1 ). (2.9) 

In terms of the function W definition of uniform material (|2.8|) takes very simple 
form 

W{X,A) = W(Y,A). (2.10) 
Thus, the uniformity condition (|2.4[) for the strain energy function VT is 
equivalent to the statement that the function W(X, A)), X e M,Ae GL(V) 
does not depend on the point X 6 M. As a result, W(X, A) it is a function on 
the linear group GL(V) only. This result is the central point of the theory of (first 
grade) uniform hyperelastic materials. It reduces the study of material properties 
of body M and the evolution of those to the study of the uniformity mapping Px 
and the function W on the linear group GL(V). 

Additional physical requirements (e.g. material frame indifference, presence of a 
nontrivial material symmetry group, etc.) lead to additional restrictions on the form 
of the energy function W. For instance, material frame indifference requirement 
leads to the conclusion that VK(A) is a function of invariants of matrix A. If a 
uniform material is isotropic, function VK(^4) is left invariant with respect to the 
multiplication by elements of 50(3) ( [40 ] ITT ] [34] ) . 

Returning to the the strain energy density function W(X, F(X)) we see that for 
a uniform material with the uniformity mapping P the strain function W takes the 
form 

W(X, F(X)) = J P {X)W{F(X) o P{X)). (2.11) 

2.6. Material metric of a uniform structure. As it was already known to E. 
Cartan (see [3]), to a zero curvature linear connection u> (absolute parallelism) 
determined by a frame p (or by the corresponding uniformity map P) there cor- 
responds the material metric g defined as the pullback of Euclidian metric h by 
the mapping P^ 1 

g(A) = P x lh. (2.12) 

This definition is equivalent to the declaring the frame p g-orthonormal at each 
point X s M . In local coordinates X 1 the metric g has the form 

giJ {X) = (Px'TAPx 1 )^ = {D{X)^)f (D(X)-^g MN , 

the first expression being given in terms of the uniformity mapping P while the 
second is in terms of the corresponding (l,l)-tensor field D. 

The curvature of the metric g is then defined by the torsion of the connection to 
(see [221 [39]). 

2.7. Examples. Elastic strain tensor of a body in a configuration <f> is defined by 

E? = i/n(g - 1 C(^)) « ^(g^CW - I), 

where second expression is the linear approximation of the first one, 021 E2]- Recall 
that the strain energy function of an isotropic material in linear elasticity has the 
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form 

W{4>) = A[Tr(E e ')] 2 + /xTr[(E eZ 2 )], 
where A, \i are Lame coefficients (see [39]). 

Using the same function W — \[Tr(A)) 2 + [iTr[A 2 ] on the linear group GL{V) 
but a nontrivial uniformity mapping P, we come to the model of a quasi-isotropic 
material. Uniformity mapping P defined the material metric g as above. This 
allows to redefine the elastic strain tensor using metric g instead of the reference 
metric go : 

E ei = ~ln(g- l C(0)) » \{g^C^) - I). (2.13) 
Strain energy of a quasi-isotropic material in linear elasticity is defied as follows 

W P (X, F(X)) = » P (X)[\(Tr(E el )) 2 + ^Tr{W l 2 )]. (2.14) 
It is easy to see that the strain energy is the quadratic function of the conventional 
elastic strain tensor with the tensor of elastic moduli depending on material 
point X. 

Another example is provided by a quasi-Hookean material (see [22], p. 11), i.e. 
the uniform analog of the neo-Hookean material with 

W{4>)=a[Tr(E* 2 ) -3]. 

The quasi-Hookean material corresponding to a uniformity structure P is defined 
by the same strain energy function but with the redefined strain tensor E e ' = 
7}ln(g~ 1 C((j>)) where material metric g is used instead of the reference metric go 

W P (X, F(X)) = a(Tr(E ei 2 ) - 3)]. (2.15) 

In the case of a homogeneous uniformity structure last expression reduces to the 
strain energy of standard neo-Hookean material. 

2.8. Evolution of the uniform structure. Evolution of the properties of a uni- 
form material is characterized by the time-dependence of the uniformity mapping 
P and that of the function W. An appropriate characteristic of the evolution of 
uniform structure P is the material velocity L(X) that was studied by different 
authors, see for instance [29l[12l [T]. 

The material velocity of the uniformity structure P is defined as the material 
point and time dependent linear mapping 

L t (X) = P x 1 o^:V^V. 

Under a loading both the uniform structure P and the deformation mapping <fi 
are evolving. As a result the couple (Px(t), 4>(t, X)) (or (g(t, X), <j)(t, X)) describes 
both the (total) deformation of a material and the evolution of its properties (elastic 
moduli, reference density, etc.). The rate of change of this couple is given by 
(L t (X),V(t,X)), where V(t,X) = is the physical velocity. 

3. ELASTO-PLASTIC MULTIPLICATIVE DECOMPOSITIONS OF THE DEFORMATION 

GRADIENT 

At the end of 1950s B.Bilby, E. Kroner ([17] ) and later on E. Lee ([18]) pro- 
posed the following multiplicative decomposition of the deformation gradi- 
ent (BKL-decomposition) 

F = F e F p (3.1) 



8 



V. CIANCIO, M. DOLFIN, M. FRANC AVIGLIA , AND S. PRESTON 



as the product of two smooth (l,l)-tensor fields of elastic and plastic defor- 
mations, respectively. To provide a geometrical illustration of this decomposition 
an intermediate configuration C t * was introduced between the material body 
M and the current configuration Ct — 4>t{M). 

Decomposition F = F e F p is used to study the behavior exemplified by an elasto- 
plastic behavior of a material which undergoes deformation under a slowly applied 
load beyond the elastic range and then, after unloading, preserves some "perma- 
nent" strain (deformation) . We refer the the monograph [53] for more examples and 
references concerning multiplicative decompositions of the deformation gradient F 
and their applications. 

3.1. Relation between the BKL-decomposition and the theory of uniform 
materials. Recall that the deformation gradient F t (X) of a configuration : M — > 
E 3 is the two-point (l,l)-tensor field in M defined by the linear isomorphism of 
the tangent spaces 0») : T(M) -> T((f> t (M)) at X G M. Here C t = cf> t (M) is the 
configuration of the body at the time t. 

The decomposition (3.1) can be hardly interpreted other then as the composition 
of tangent bundle mappings over some mappings of corresponding base manifolds 

T(M) ► T{C* t ) ► T{cj> t {M)) 



M ► C* t ► <j>t{M) 

since the tensor fields F e and F p should be strictly anchored at some manifolds 
(domain and target of each). Moreover the first mapping F p should define a map- 
ping from the tangent space Tx (M) at a point X £ M to the tangent space at some 
point Yx of intermediate configuration C* . The correspondence X — > Yx should be 
one-to one, otherwise the composition (3.1) cannot be an isomorphism of the tan- 
gent bundles. Therefore, there exists a unique one-to one mapping 0i : M — > C\ 
underlying the tangent bundle mapping F p . Mapping 0i can be assumed to be 
diffcrcntiable. 

In the same way F e can be viewed as a mapping of tangent bundles T(C ( *) — > 
T(Ct) over the differentiable mapping 02 : G\ —> T(Ct) of basis manifolds. 

We obviously have <fi — <p2°4>i- Therefore, 02 is onto. Restricting, if necessary the 
intermediate configuration manifold one may assume, without loosing of generality, 
that 0i is onto and <p2 is one to one. Thus, both 0i and <p2 can be considered as 
diff eomorphisms . 

Remark 3. Defining the decomposition (3.1) some authors presume that the map- 
pings F e and F p are nonsmooth or even noncontinuous, reflecting microdefects 
densities in the manifold M. Translating this into the language of tensor fields 
and using the derivatives of these tensor fields one should however assume some 
smoothness. Usually it is done by considering these tensor fields as smooth averaged 
characteristics of the structural state of the material. 

Remark 4. Mapping <f>\ presents the intermediate configuration introduced in 60th 
by a variety of researchers, see [TH |37j [35] . It was used for the construction of plastic 
deformation gradient F p and the elasto-plastic decompositions of total deformation 
gradient F but as far as we know, was not considered previously as an independent 
dynamical variable. 
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Now we are ready to make the next step. Consider the tangent mapping 0it* : 
T(M) -> T(C*) and compare it with the mapping F p (t,X) : T(M) -> T(C t *). 
Since mapping 0i t * is linear isomorphism at each point X € M , one can write, for 
all tangent vectors £ £ Tx(A^), 

F*(t;(*,O)=0w.x-oD t (JO.£ 

where D t (A) is uniquely denned smooth (l,l)-tensor field in M. 

In exactly the same way one can present 

F e (t,y,r;)=<WoF«(<,y).£ 

for the uniquely defined smooth (l,l)-tensor field F e *(t,Y) in C t *. 

If we pull back the (spacial index of) tensor held F e * (t, Y) from C t * onto M by 
the differential 0i t * of the point mapping 0i we get another (l,l)-tensor D e on M. 

Since <j>{t,X) = t2 °0ti for all t and, therefore, (j>(t,X) m x = 4>t2 *^ tl (x) °0ti *x, 
combining this with the decomposition (3.1) we get 

^, = F e oF p = (0 2 , o D e ) o (0!, o D) = (0 2 * o las ) o o D e o o D = 

= (02 o 0i)„ o (D e o D) = 0, o (D e o D), 

so that 

r> e (t 1 x)-r>(t 1 x) = id T{M) . 

As a result, being transferred to the material manifold M, the (l,l)-tensor 
helds connecting integrable mappings (j)^ (i — 1, 2) to the tangent bundles mappings 
F p and F e are inverse to one other. This is hardly a surprise since in the physical 
literature only one of these tensors was considered as an independent dynamical 
variable; see [6l [26], 

In the same way, 02 = ° <\\ x would be also redundant. As a result, the only 
independent dynamical variables in this scheme are diffeomorphic embeddings 0, 0i 
and the material (l,l)-tensor held D. 

Remark 5. One can of course choose another triple of variables as independent 
dynamical quantities, for instance one may use (0i, 02, D) if it is preferable to deal 
with the elastic deformation 02 explicitly. 

Remark 6. We consider here only the decomposition F e F p , but the same argu- 
ments would produce a geometrical representation of the reverse F p F e -decomposition 
as well. 

Remark 7. Notice that the choice of an intermediate configuration (C t *, 0i t) par- 
ticipating in the decomposition (|3.1jl is far from being unique. In particular, let 
us show that we may formally choose the image C t * = O (M) as the intermedi- 
ate configuration with 0i t = O being time independent. To do this denote by 
tfj : C t * — + C the diffeomorphism ip — O o 0^ 1 . Transfer the tensor F e to C a 
as follows: F e = "0»(F e *), where we are using the diffeomorphism ip together 
with its inverse to push forward the (l,l)-tensor F e *. Thus we get the mapping of 
tangent bundles 

FF C = X t*°F eo = X t* °ipt*oF e *o 4>~\ 



10 



V. CIANCIO, M. DOLFIN, M. FRANCAVIGLIA, AND S. PRESTON 



Define also the diffcomorphism \t ■ C Q — > Ct as \t — 4> fa. 1 • an d define the 
mapping of the tangent bundles by setting: 

F p = OT oF pM . 

Then we have (j> — x ° fa, and, as is easy to check, F e o F p = fa as required. 

As a result we get a simplified scheme of the elastic-plastic F e F p -decomposition 
of the deformation gradient F = fa of a uniform material. Notice that the in- 
tegrable part fa t of the plastic deformation gradient F p is lost in this simplified 
scheme. That is why it is preferable to work with the previous scheme where the 
intermediate configuration is different from the image of the reference embedding 

fa- 

Remark 8. Notice that the couple (D,fa) represents another model of evolution 
of the material of the same type with the same uniformity structure. This 
model of pure inelastic evolution is related to the model (<j>, fa , D) by the elastic 
deformation fa. 

If we start with a time dependent uniformity mapping P t and two configurations 
fa fa : M — > R 3 , then one can (reversing the arguments above) construct the 
"elastic deformation" fa — 4> o fa 1 and the mappings of tangent bundles F p : 
T(M) -> T(C* = Im(fa <t ), F e : T(C* = Im{fa, t ) -> T(C t = Im(fa)), such that 
the construction above returns us to the triple (Pt,4>, fa). 

Finally, there is a freedom in the choice of the decomposition F p = fa* o D given 
by an arbitrary diffcomorphism ip € Diff(M): 

F p = ^oD=(^ 1 o^ 1 ),o(^oD). (3.2) 

Thus, we can introduce the following equivalence relation between the pairs 
(fa, D) of the (time dependent) mappings fa : M — ► R n and nondegenerate (1,1)- 
tensor fields D in M. We say that two pairs (fa, D), (xi,K) are equivalent if there 
is a diffeomorphism tp G Diff(M) such that 

X i=fa°ip-\K = ^oV. (3.3) 

Collecting the considerations presented in this section we get to the following 
conclusions 

(1) BKL-decomposition F — F p F e of the deformation gradient F = fa, of a 
(total) configuration fa : M — > E 3 presupposes the existence of (inter- 
mediate) inelastic configuration fa t ■ M — > E 3 and of the non-degenerate 
(l,l)-tensor field D t in the material space M such that 

F p = fa, o D, (3.4) 



F^^oD-'ofc 1 . (3.5) 

(2) Configuration fa is the mapping fa t : M — > E 3 defining the integrable 
part of inelastic (plastic!) deformation gradient F p , 

(3) (l,l)-tensor D is equal to 
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(4) Vice versa, to any triple (fa, Ff,Ff) consisting of two configurations and 
a non-degenerate (l,l)-tensor field Dt in M there corresponds the multi- 
plicative decomposition 

fa = F = F e o F p 

of the deformation gradient F of the total deformation history fa 
Tensor field D t defines the time-dependent uniform structure P t = D t o P 
(see Proposition l)in the material body M. 

Uniform structure D determines the time dependent Riemannian metric g t 
in the material manifold M (see Sec. 2. 6) by the formula: 

gt u = hijOiDfj. (3.6) 

Remark 9. Notice that the decomposition F = fa = F e oF p determines the tensor 
field D and the plastic integrable deformation fa up to an action of a diffeomor- 
phism x- 

On the other hand there are arguments showing that the inelastic configuration 
fa is defined uniquely by the history of deformation. If total deformation fa — > fa 
is going under certain conditions of loading, heating etc, the unloading or turning 
of the heat at some moment t\ produces certain configuration fa : M — > E 3 . Often 
the unloading happens fast and is not accompanied by an essential change of the 
material structure. See [19) where relaxation of a material to the intermediate 
configuration fa during unloading is discussed. As a result we may associate with 
the moment t\ the final configuration fa tl taken by the body M after unloading. 
As a result, under physically reasonable assumptions on the evolution process, the 
intermediate configuration fa, and therefore the tensor field D, are determined 
uniquely (up to a composition of fa with the Euclidean motion of the physical 
space E 3 ). 

Using the terminology of G.Maugin and W. Muschik ([28 ), the configuration fa 
is "observable but not controllable'''' (or, citing the same work, "partly controllable" 
through the loading conditions). 

In this analysis of the kinematics of the BKL-decomposition we identified three 
variables: the uniform structure P t , the integrable part of the "plastic defor- 
mation gradient" fa (intermediate configuration), and the total deformation cj), 
which we will consider as physically independent. Thus, the full dynami- 
cal/thermodynamical picture with this kinematics should include all three com- 
ponents. Below, discussing the dynamical structure of the presented model, we will 
be using material metric g instead of the tensor field D as the material dynamical 
variable. Metric g contains the essential information about the uniformity struc- 
ture P t and is more convenient to use when describing the evolution of material 
properties then the tensor D or uniformity mapping P, see |11) . 

4. 4-Metric Theory and Strain Tensors 

4.1. Strain Tensors. Nonlinear Elasticity Theory has, as its geometrical keystone, 
the question of the comparison of two metrics: the reference metric g G and the 
Cauchy metric C(fa) of a configuration fa, [22] . 

As Theorem 1 shows, the multiplicative decomposition leads, in its geometrical 
form, to the presence of four metrics in the material manifold M: the reference 
metric g Q , the material metric g generated by the (l,l)-tensor field D = F p M , the 



12 



V. CIANCIO, M. DOLFIN, M. FRANC AVIGLIA , AND S. PRESTON 



Cauchy metric of the integrable part of plastic deformation C(<f>i) and finally, the 
Cauchy metric of the total deformation C (</>). It seems natural to define appropriate 
strain tensors as measures of comparison between pairs of these metrics and use 
these tensors for description of different processes developing in the material. 

We will introduce six strain tensors as suited to describe the state of our solid 
and characterize specific processes undergoing in the body. 

(1) Elastic strain tensor 

E e < = l -ln[C{^)-'C{4>)\ « l -C(^[C{4>) C(&)] « Iz'm) C(0O]. 

Elastic strain tensor measures the elastic part of the deformation at each 
instant of time and vanishes under unloading. Tensor E el and its linearized 
variant are C(</>i -symmetrical (i.e. Efj — C{<px) IK E e j K is symmetrical 
(0,2)-tensor. 

(2) Inelastic strain tensor 

E m = i?n[g- 1 C(^ 1 )] « \z\C{4>) - g] « \&\C{4>) g]. 

Inelastic (plastic) strain tensor measures the plastic but still Euclidean 
deformation of the body i.e. permanent after unloading but not leading 
to any residual stresses in the material. Tensor E m and its linearized 
version E™ n are (/-symmetrical, i.e. (0,2)-tensors gixE™ K ,giKE™ n K j are 
symmetrical. 

(3) Material strain tensor 

E ro = ^n(gJ 1 g)«^ 1 (K-8o)- 

Material strain tensor measures the pure metrical evolution of the material, 
not leading to any deformation (material points displacement). Tensor E m 
and its linearized version E™ n are go-symmetrical. 

(4) Euclidian strain tensor 

E eud -i/nfe-'CWl^^g-MCW-g]. 

Euclidian strain tensor measures the integrable part of the total deforma- 
tion. 

(5) Total strain tensor 

E* ot - ^n(g" 1 C(0)) « ig- 1 (C(0) - go ). 

Total strain tensor measures the decline of the Cauchy metric of total defor- 
mation <f> from the reference (euclidian) metric g Q . It is observable. Tensor 
E tot and its linearized version E*°£ are go-symmetrical. 

(6) Total inelastic strain tensor 

E* m = ^^(g-^^)) « ^g-^C^r) - go ). 

Total inelastic strain tensor measures the decline of the Cauchy metric of 
inelastic deformation tfii from the reference (euclidian) metric g Q . It is 
observable (after unload). 
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In each case we provide the linear approximation form(s) of the strain tensors suited 
for small deviation of the former metric from the latter one. 

Remark 10. In some simple cases, say when the (l,l)-tensors g _1 C((/>i) and 
C(0i)- 1 C(0) commute, from the relation g _1 C((/>) = g^C^i) • C(0i)- 1 C(0) 
we conclude that the linearized Euclidean strain tensor E etici splits as follows 

E eud „ lg-l[ C (0) - g] = ig-^Cfa) - C(0x)] + ig-MC^l) - g] « 

w g- 1 C(</-i)E e/ + E m w E e/ + E™. (4.1) 

If all three (l,l)-tensors g ( 7 1 g, g _1 C((/>i) and C(0i) _1 C(</>) commute we have 
an additive decomposition of the linearized total strain tensor E tot 

E tot » igo '[CW - go ] = -g^TO) - C(^) + C(0!) - g + g - go] = 

= ig o - 1 C(0 1 )C(0 1 )- 1 [C(<^) - C(^)] + igQ-^g-MC^i) - g] + ig^fe - go) = 

= E m + g^gE™ + g„ x C(0i)E e! « E m + E m + E ei . (4.2) 

4.2. Choice of the dynamical variables. In order to determine which combina- 
tions of dynamical variables (</>, <f>i, D) and their derivatives might enter the internal 
or free energy, dissipative potential, entropy and other dynamical and thermody- 
namical quantities we have to take into account requirements of invariance or co- 
variance of these quantities with respect to the appropriate material and spacial 
transformations. For instance, the frame indifference requirement (1321 131] ) leads 
to the conclusion that the deformation gradient F(X) of the total deformation <f> 
enters these quantities only in combinations C(0i) -1 C (<j>), g^ 1 C{(j)) or g~ 1 C(cj)). 

Material metric g and the Cauchy-Green tensor C((f>i) depend on the choice of 
the plastic decomposition (|3.4[) . Thus, it is important to determine which tensors 
or combinations of tensors constructed from the dynamical fields (4>, 0i,D) are 
independent of the choice of the plastic decomposition (|3.4p . 

We consider several such combinations. 

1. The total strain tensor 

E tot = iln(gjV*h) « ig- 1 [C(0) - go] (4.3) 

is independent of the plastic decomposition (|3.4|) . Since it is a (l,l)-tensor its invari- 
ants are independent of the decomposition (13. 4[) and, moreover, one can combine it 
with other similar tensors to produce new invariant combinations. 

2. The plastic deformation gradient F p = cf>i* o D is invariant under decomposi- 
tion (|3.4p but it is convenient to transform it into a material tensor. For instance, 
one can use the following variant of the Cauchy-Green tensor 

F p *h = D*(<^h) = D*C(0i). (4.4) 

Lifting one index in this tensor by means of the reference metric go we get the 
material (l,l)-tensor g" 1 *F p *h. This tensor carries information on both integrable 
and nonintegrable parts of the plastic deformations. 

3. The material metric g = P _1 *h can be written in a number of different ways. 
For instance, by using the reference metric g D = P _1 *h we can get 

g = p-^h = P- 1 *(P *g ) = (P o P- 1 )*g = (P o P- 1 )- 1 * go = D- 1 * go . (4.5) 
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In local coordinates we have 

— r n-i\M r r>-i\N „ „-i ab _ n A r> B „ MN (a r\ 

QAB - \J-> ) A \ U )b9o MN, 9 - Um U N% ■ ( 4 -0) 

Under the change of decomposition (|3.4p , the material metric g is transformed 
by a diffeomorphism i/j into a new metric g 1 ^ 1 as follows 

g^ = o D)- 1 *g - (D- 1 o ^)*g - i?- 1 * o D 1 *g = rs- (4.7) 

In coordinates 

{ji*)ab = {r x )%{r x )%9us, u* -T D = tfri&g- 1 MN - 

Therefore, the curvature tensor Rj KL (as well as the corresponding Ricci tensor 
Rjj(g)) is transformed tensorially by ip* or ip*. Thus, invariants of these curvature 
tensors (in particular the scalar curvature i?(g)), do not depend on a choice of 
plastic decomposition (|3.4|) . 

4. At the same time the Cauchy metric of the inelastic deformation C(</>i) is 
transformed by tp -1 tensorially as well: 

Cfa of 1 )^ = (ip- 1 *o,pih) AB = [V" 1 *C(0i)]ab = {^a^'^bC^mn. 

(4-8) 

Combining the material metric and the Cauchy tensor of the inelastic deforma- 
tion (f>i we finally get the inelastic strain tensor E m . For a different decomposition 
(|3.4p we have 

= ^g- 1 ^s^^-YsCiMKL = KA^Bg- 1 mk c(Mkl = 

= <(^- 1 )l(^ 1 C(0i))f. (4.9) 

Thus, E m transforms tensorially under a change of plastic decomposition (|3.4p 
and its invariants are independent of this decomposition. 

The presence of the two quantities F p and E' m whose invariants are independent 
of the choice of plastic decomposition makes it important to compare these two 
quantities. We have 

g-^F^h) = g - 1 (D*C(^ 1 )) = g; AC D%D2C(<h)MN = 

A (g-' KC Dk Dq )C (4>i) m n D B = A (g^ ML )C(<f>i) MN D» = 

= Dl 1 A (g^C{<l> l )) N D N B = Dl 1 A (exp(2E m )) L N D N B , 

i.e. 

g-^F^h) = D~ L l A [exp(2E m )] L N D^, (4.10) 

and therefore invariants of the tensor g ( 7 1 (F p *h) (i.e. tensorial characteristics of 
the plastic gradient deformation F p ) contain the same information as the invariants 
of E in . 

Thus in the restricted case it seems natural to choose an internal energy it as a 
function of the invariants of the two strain tensors E to * and E m , of the invariants 
of the curvature tensor Rj kl (g) (of the Ricci Tensor i?ic(g) j in the 3-dim case, see 
01), temperature and its g -gradient: u = u([E tot (0), E™(0i, g), Ric(g), 9, V s °6»]). 

If we adopt the assumptions of Remark 5 (i.e. removing the restriction to use 
only tensors invariant under the plastic decomposition (3.6)) we may consider all 
three strain tensors E tot , E m , E m (one can replace E tot in this list by E e * if it is 
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preferable) as independent dynamical variables and, together with the Ricci tensor 
of the material metric g, include them as arguments in the internal energy: 

u = u[E el ,F m ,F m ,Ric(g),9, V s °(9]. (4.11) 

In this approach the effects of different types of processes are directly separated. 

4.3. Additional strain decompositions. Between the strain tensors introduced 
above, conventional strain tensors and both deformation gradients F e , F p of the 
multiplicative decomposition there exist different relations that may be in some 
partial cases more convenient than those presented above. Below are two examples 
of such relations, the first being valid in the linear case, the second in a nonlinear 
situation. 

1. Linear case: 

Since F — 0* = F e o F p we have 

g"V*h = g^F'fF'h - h + h) = g-^h + g- x Ff*(F e *h - h) 
and therefore 

Kn = i(g-V*h-S) = ifej^h-tfJ+g^F^F^h-h) « D^E&D+g^I** 

(4.12) 

where we defined 

1 

2 ( 

This decomposition can be compared with those in Section 4.1. 

2. Nonlinear case: 



E^ n = -(F e *h-h). 



E tot = ^M^fh) - ^nCg-^F^hjh-^F-h)) = i^(g- 1 (F^h).exp(2El)) 

^(g-^F^h) ■ (F?*exp(2Ef ld ))) « ^nfeo ^h) + g^F^E^), (4.13) 
here F e J ld = ±ln(h.- 1 F el *h) as in conventional finite elasticity. 

4.4. Strain Rate Tensors. We define the strain rate tensors as time derivatives 
of strain tensors. As a result we get the strain rate tensors E tot ,E m ,E e ' and E" 1 . 

On the other hand there are other rate characteristics for each of the three 
participating structures, i.e.: 
(1) The material velocity: 

h D (t,X) = B(t,X)- 1 o^D(t,X). 

This velocity is related to the velocity Lp = P^ 1 o P t introduced in Sec. 2. 8 
by the relation Lp = P o Lp o P^ 1 . We see now that D = DLp and 
therefore we get the relation between the material strain rate tensor and 
the material velocity Lp, used in (TJIHIHH] and other papers. 

E m = \%- l k = ig- mn + D?D$g MN ) = 

— 2 So \ U K L D l U J 9o MN + L> 1 U s L D jg MN ) - 

= l^o 1 (Ld i9p kj + L S D jgis) = g-^Ln)*! ■ g mj) , (4.14) 
where in the last formula there is symmetrization by indices IJ. 
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(2) The total velocity is defined as 

v(x,t) = ^4>(x,t) 

and its gradient is related to the linearized total strain rate tensor E'°^ = 
ig o " 1 (C(0) - g Q ) by the relation 

(ffoEj^)MJV = hijV^ M <t? N y 

Being written in Euler (spacial) coordinates this relation reduces to the 
standard one ([55]). 

(3) Finally, the velocity of the inelastic deformation 

is related, in a linear approximation, to the (linearized) inelastic strain rate 
tensor ~E\f n by a relation containing the symmetrized velocity gradient and 
the material velocity L^. In the calculation that follows we are using the 
formula JjC(</>i) = 2h i jV^ M 4> J 1 ^ for time derivative of Cauchy-Green 
tensor of configuration <f>\ . We have 

Et = |(ig- 1 W 1 )- g ))=^(g- 1 -C(0i)) = 

= ig-i^cwo + ^(g- 1 ) • cc^) = g-M^ntM^U)] + ^g- 1 ^ 1 • c(0t) = 

= S-HhijV^M^} +g- 1 (g o E ro )g- 1 C(0 1 ) « (V g ° Vl ) sym + E m . (4.15) 

Here Vi A ,/ = ^M^-ij^i is the covariant form of the convective velocity of 
inelastic configuration </>i , see [35] . In the last approximation we replaced 
g ~ g in the first and the second terms and C((f>i) ~ g in the second term. 

5. Lagrangian, free energy, dissipative potential and the stress 

TENSORS 

Dynamical equations describing the evolution of the system characterized by the 
variables (<fi, 4>i > g) are obtained by combining the canonical (Lagrangian) compo- 
nent and the the dissipative forces. 

The Lagrangian in our model is the combination of kinetic, potential and internal 
energy terms: 

L = K- Pref u - U{<j>), (5.1) 

where K = | |V|^ is the density of kinetic energy, U{4>) is the potential of the 
volume forces and u is the internal energy per unit of mass (see 4.2). 

5.1. Lagrangian and internal energy. It is traditional to define free energy 
density -0 as a function of the elastic deformation gradient, the temperature a 
material point X and additional internal parameters a (see [24j[28], Sect. 10. 1A). 

Dissipative pseudo-potential is, in this approach, the function of rates of de- 
formation gradients and time derivatives of internal variables T> = D(F e , F p , a, Vi?, Va, $) 
([24l [28]). This allows one to define the total, elastic and plastic stress tensors and 
the thermodynamical forces conjugate to the parameters a, thus separating differ- 
ent factors in the dissipation inequality (see 10.21 in [25], Sec. 10.1 or [24]). 
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Comparing the expression for the internal energy (|4.11j) with these in [211 [28] 
one sees that the metric g entering the free energy through the tensor E m plays 
here the role of an internal variable a and its Ricci tensor Ric(g) takes the place of 
the space gradient Va ([HI [28]). The elastic E e/ (respectively inelastic E m ) strain 
tensors are direct material analogs of e e (respectively of e p ). Thus it is conceivable 
to adopt the internal variables approach in searching for the form of the equations 
governing the behavior of our system. 

We take the Lagrangian in the form 

L = L(p ref ,V, E e <, E m , E"\ ittc(g),0, V 9 °ti) =K- p ref u - U(<j>) = 

= ^|v| 2 ,+p re/ [7.ffi c (g)+^||^ 

(5.2) 

In this expression 7 is a constitutive tensor, U is the potential of body forces, /o 
is the "basic inelastic energy" and / is the strain energy of linear thermoelasticity, 
i.e. a quadratic function of the elastic strain tensor with coefficients depending on 
the temperature and the remaining inelastic strains: 

/(E in , E m , 1?; E ei ) = [c + (0 - 0o)ci] : E et + (e : E eZ : E ei ). (5.3) 
Here Co and Ci are tensors characterizing the interaction of the elastic processes 

with the inelastic ones and temperature respectively (for instance, Ci is the thermal 

expansion tensor). The tensor e is the elasticity tensor. 

Assuming that the decomposition (4.2) is valid, substitution of the total strain 

tensor E tot instead of the elastic strain tensor E e/ into the expressions for internal 

energy and dissipative potential (below) does not change the form of function (5.3) 

/(E in , E m , tf; E ei ) = /(E in , E m , tf; E tot - W n - E m ) = 

= [c o +(0-i?o)cx] : (E*°*-E in -E m ))+(e : (E*°*-E in -E m ) : (E to *-E" l -E m )) = 

= +(e : (E™ + E m ) : (E™ + E m )) + 

[(c - 2e : (E™ + E m ) + (1? - i? )ci] : E* ot + e : E* ot : E tot , (5.4) 

but changes the tensor Co and adds a term to the inelastic energy Jo- 

This allows us to replace E e ' by E tot in the internal energy, so that we can 
equivalently use Lagrangian in the form 

L = L[p ref , V, E tot , E m , E m , Ric(s), i5>, V 9 °^] = 

(5.5) 

The strain energy here has the form (5.2) where E ei is replaced by E tot . 

In the dynamic equations it is more convenient to use E* ot since <fr is a geometri- 
cally explicit and observable quantity while in the dissipative inequality it is more 
convenient to use E ei because it allows one to separate inputs of different processes 
into the entropy production. 

5.2. Free energy. Free energy is defined as usual, by the equality 

ip = u — s6, 

where s is the specific entropy. 
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5.3. Dissipative potential. The dissipative (pseudo) potential id chosen to 
be a function of the following variables 



V = £>(E"\E m ;E m ,t?). 



(5.6) 



We include E m together with g to emphasize the difference between the kinetic 
energy related with 0* ot and the inelastic strain rate E m participating in the process 
of viscous dissipation, see [27] . 

5.4. Evolution equations. Introducing the action for the material in a domain 
G C M , corresponding to the Lagrangian L 



Ldv n 



G 



we will write down the equations of evolution for the system characterized by the 
dynamical variables <f>, (f>\, g: 

(1) Equilibrium equation: 

S i- <") 

Since the total deformation <f> enters only through the elastic strain tensor 
E ei this equation is, essentially, the equilibrium equation of elasticity theory. 
If the strain energy is chosen in the form (5.3), this equation takes the 
conventional form of the elasticity equilibrium (Euler) equation with the 
elastic moduli depending on the material point and temperature (see (5.3)) 



1 d(p re f\/\g\hi 3 V :i ) 
dt 



P 



(5.8) 



Here u = —dU (cf>) is the 1-form of the body forces and P/ is the elastic 
first Piola-Kirchhoff stress tensor (see below Sec. 5. 5). Covariant derivative 
is taken with respect to the material metric g. 
(2) Equation of plastic deformation: 

SA SB , , 

m = ~w (5 ' 9) 

Notice that in difference to the usual form of this equation (\24\ I28j) we 
take the variation of dissipative potential by the derivative <f>\ of internal 
variable <p\ rather then the partial derivative. This is necessary due to the 
fact that <pi enters Lagrangian through its spacial gradient. 
To clarify the form of this equation we notice that 



SA 



1 



du 



-div g {V{ 



i) l.-lh M 
rriin N U£j N 



1 M 



'1,1 



T]S N 9 MK h ij W N S I K + ^S I N ) . (5.10) 



Here V[ 4 is the first Piola-Kirchhoff stress tensor density of inelastic con- 
figuration 4>i, see below. Last equality is due to the fact that N — 
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On the other hand, 



SD 



= -d 



K 



an dEp M 



8E™ M 



l.K 



= -8 



K 



8T> 



8E$ M 



\ h ij{9 



MKii 



■9 MS 4,s^) 



where we have used (4.15). As a result, equation (5.9) has the form 
<9D I 



diVgiPij+d. 



K 



8E™ M 



1,N 



■9 



0. 



(5.11) 



(3) Equation of metric evolution: 

5 A _ 8T> 
Sg 8g ' 

Here we have used the partial derivatives in the right side of equation since 
i/n(g " 1 g) depends on g but not on its derivatives. 



(5.12) 



E'' 



If the free energy depends on the scalar curvature R(g) instead the full 
Ricci tensor, equation (|5. 12|) has the form 



c-£(gY J =-T' J - 



8T> 

<9E m 



(5.13) 



where 



T 



I J 



([^|V| 2 , + PT*f\lAV*'i% + fo + /] - CTWlv^gl) 



V\s\ S 9u 

is the canonical energy-momentum tensor including elastic effects, effects of 
inelastic deformation and some thermal effects and £{g) IJ is the Einstein 
tensor of metric g ( [311] ) • If g = go is the reference metric then this equation 
is absent (go is fixed). 



Remark 11. 



In the 2-dim elasticity any metric g in M is Einstein metric, i.e 
In this case using the scalar curvature R(g) instead of the Ricci 
tensor in (5.13) does not place any restrictions to the material metric g. 



5.5. Stress tensors. Stress tensors characterizing material response to the defor- 
mations, heating and other physical processes play a crucial role in the formulation 
of the evolution equations and dissipative inequalities. In the material (Lagrangian) 
formulation there are several stress tensors playing different roles in the dynamical 
picture. They are related to one another and, through the deformation 0, to the 
only stress tensor that is usually present in the Euler picture - the Cauchy tensor 
<t. Such a plurality of material stress tensors is related to the presence of two ma- 
terial metrics, i.e. g D , g, used to raise and lower the indices in tensors and the two 
different Cauchy metrics - C(0), C(4>i). 

Here we recall the definitions of the most useful stress tensors through the total 
internal energy u or the strain energy / and provide the formulas relating them to 
one another. 

For the total deformation ip we introduce three stress tensors defined by the 
differentiation of internal energy by the deformation gradient F = Cauchy 
metric C{(j>) and the strain tensor E tot . 

The formula relating S and T el is obtained in the assumption of linear approx- 
imation C(4>) — C(4>i) + 2C(cf>i)E el , see Sec. 4.1. J(4>) here is the Jacobian of 
the total deformation <fi calculated with respect to the metrics h and g , see |22j . 
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Type 


I Piola-Kirchhoff 


II Piola-Kirchhoff 


Strain dual 


Tensor 


•pi — n f 9u 
I i — Pref 


S — tyref qc{$)ij 


Tel I _ „ . du 
1 J Pref QEf J 


Relations 




v H = S IK ^ K 


Tf 1 = S^Ci^Kj 



Table 1. Stress tensors defined by total deformation 4>. 



Sec. 2. 2. Expression for the tensor S tot IJ in the Table 1 is the material Doyle- 
Erickson formula (see [23]). 

For the inelastic deformation cf>i we introduce three stress tensors defined by the 
differentiation of internal energy by the deformation gradient F = </>i», Cauchy 
metric C((j>i) and the strain tensor E tm (see Sec. 4. 1): 



Type 


I Piola-Kirchhoff 


II Piola-Kirchhoff 


Strain dual 


Tensor 


n i =Pref^- i 


S\ %Pref 8C(0)i j j 


rpin I _ du 
1 J Pref gE inJ 


Relations 




= Sl K <f>\,K 


^1 9kj-1j 



Table 2. Stress tensors defined by inelastic deformation <f>\. 



The relation between Si and T m is obtained in the assumption of linear approx- 
imation E m = icrHcOi) - s)- 

For the deformation (evolution) of material metric go — > gt there are defined the 
following stress tensors 



Type 


Eshelby stress, QH EE1 [25] 


Canonical 


Strain dual 


Tensor 


hi — n 9u 

— PrefgpT 
i 




T*m I „ du 


Relations 




b) = S m MN P- X l g NI 


& go kj — 1 j 



Table 3. Stress tensors defined by deformation of the material metric. 



Here f is the uniformity mapping P : V — > Tx{X) and the internal energy 
is refereed to the reference volume d go V, [IT], Ch. 5, Sec. 5. 5. Formula relating S m 
and T m is obtained in the assumption of linear approximation E m — \g Q ~ 1 (g — go)- 

One can define the variant of the Eshelby stress by bj = —po-§§r using the 

material (l,I)-tensor D. Its relation to the tensor b\ is given by bj — b l jP ~l 1 . 

Notice also that the canonical stress tensor S m is the direct analog of the spacial 
part of the energy-momentum tensor of the Gravity Theory, [30) . 

Remark 12. It is instructive to compare our definition of the elastic first Piola- 
Kirchhoff tensor with its definition as the difference (comp. [26], Ch.10) 

Tf 1 =T[- 4 A T] n (5.14) 
where cf>2 = 4> ° <i\ X is the elastic part of total deformation. 
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6. Dissipation inequality 

In this section we present the dissipative inequality for the (cj), (fix, g)-model. 

Below V means V s ° . We will adopt here the expression (|4.11[) for the free energy 
but assume, for simplicity, the quasi-static behavior of the material (i.e. velocity 
V is negligible), potential energy U is absent and internal energy u depends on the 
scalar curvature i?(g) only, instead of on the full Ricci tensor: 

i/> = u(E el , E in , E m , iZ(g), 0, V 9 °z9) - s0 (6.1) 



with E m = 7}ln(g x g) and g playing the role of an internal parameter a (comp. 
[25]). 

We will be using the notations: 

and the formula for the time derivative of the scalar curvature 

^R(g) = lim A t^o[R(g)(t + At) ~ R(g)(t)} = £ (g) • g, (6.3) 

where £(g) = is the Einstein tensor of the material metric g (see |30|). 

Calculate now the derivative of the free energy: 



tP = T el ■ E ei + T m ■ E m + -S m • g + £(g) • g- 



V • (At?) - (V • A)0, (6.4) 



where the vectorial identity V • (Afd) — (V • A)d + A • Vi? was used and where 
tensors T in , T el ,S m are as in the Sec. 5. 

Recalling the formula for the variation — ^ — V • gyg og and using the 

notation A — wc find 

[V + 5(g)] • g = 6 / ■ g = „4(g) • g. (6.5) 
2 dg 

Then the time derivative of the free energy takes the form 

ijj = T el ■ E e ' + T m • E m + .4(g) ■ g - a ■ i? - V • (A)i9 + V • (Ai?). (6.6) 

Recall now the Gibbs inequality for a thermodynamical system with internal 
parameter a (here a = g), see |24) . 

- {ip + s0) + p t + V • (0k) - (s ■ V)0 > 0. (6.7) 

Here 

Pi = T • E* ot 

is the power of the internal work, stress tensor T will be specified below and k is 
the extra entropy flux density assumed to include contributions from the flux of 
the internal variables. 

Substituting the expression (|6.6p for ip into the Gibbs inequality (|6.7p we get 

_ T el . ge* _ T in . jw» _ J_(g) • g + ,5 . + V • (Al?) - (6.8) 

V • (Ai? - s0) + + V • (i?k) - (s • V)tf > (6.9) 
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In the special case of when one uses the linearized definitions of strain tensors (see 
Sec. 4.1) and the commutativity condition that allows to write the total strain rate 
E tot in the form (4.2) is fulfilled, the previous inequality takes the form 

(T(l + 2E m ){\ + 2E m ) - T el ) ■ E ei + (T(l + 2E m )(l + 2E el ) ~ T m ) • E m + 
+ T(l + 2E el ){l + 2E ln ) ■ E m - A(e) ■ g 

+ (s - .5 + V • A)i? + V • (0k - A9) - (s ■ V)6> > 0. (6.10) 

Now we use the fact that the derivatives E el , § are controllable variables and can 
take arbitrary positive and negative values and therefore, their coefficients should 
be equal zero, ([21]). Thus we obtain the relations 

T = T e '(l + 2E m )- 1 {\ + 2E m )-\ (6.11) 

and 

.~<,-V.A,~(g-V.*L)~g W 
Assuming for k the prescription 

k = ^ 1 A^ = ^ 1 ^, (6.13) 

the reduced dissipation inequality is obtained in the form: 

[T e '(l + 2E m )- 1 (l + 2E el ) - T m ] ■ E m + 

+ [T e ' ( 1 + 2E m ) - 1 ( 1 + 2E m )- 1 (l + 2E el ){l + 2E m )-2A(g)g ]- E m - (s ■ V)tf > 

(6.14) 

where we have used the expression E m = \(g$ 1 g — I) = ^g^g for the linearized 
metric strain tensor E m = ^g^g — I). 

Dissipation inequality (6.14) is satisfied if one request the independent fulfillment 
of the stronger conditions - two intrinsic dissipation inequalities 

f[T e '(l + 2^ m )- 1 (l + 2S e ') — T in ] -E m > 0, 

\[T el (l + 2E m )- 1 {l + 2E m )- 1 (l + 2E el )(l + 2E m ) - 2^(g)p ] ■ E m > 0. 

(6.15) 

and the thermal dissipation inequality 

-(s-V)0>O. (6.16) 

Using the relation between the tensor S m and T m presented in the Table 3 we can 
rewrite second inequality in the form 

[T el (l + 2E m )~ 1 (l + 2E m )-\l + 2E el )(l + 2E in ) - T rn - 2£{g)g ] • E m > (6.17) 

Comparing inequalities (6.15-6.17) with similar dissipative inequalities in [241 ?] 
we see that the coefficient of E m (respectively E m ) can be interpreted as the ef- 
fective stress tensor for integrable inelastic deformation (respectively, for evolution 
of the uniform structure). Such modifications of the stress tensors are customary 
in studying of the entropy production by a combination of interrelated elastic and 
inelastic processes, comp. [25], Ch.10. 
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6.1. Yield condition from dissipative inequality. If all three strain tensors in 
(6.15) are small (in comparison with the unit tensor), the inequalities (6.15) take 
the form 

trpel _ rpin] . -nun > n 

[ , J = ' . (6.18) 

[T el - T rn - 2£(g)g ] ■ E m ^ 

These inequalities can be interpreted as the yield conditions determining when 
the corresponding type of inelastic evolution (plastic integrable: <j)\ ^ and/or 
material metric gt ^ respectively) may proceed. In each case the elastic stress 
T ei should be large enough to overcome the barrier necessary for initiation of the 
corresponding process. 

This form of yield condition is similar to the condition for the plastic deformation 
to proceed obtained from the Drucker postulate, see [36], Sec. 8. 11, inequality (8.85). 

Solutions of evolution equations (5.9) and (5.11) describes also the evolution of 
stress tensors T m , T m . Therefore the conditions (6.18) for elastic stress T el evolves 
in time. This evolution can be related to the hardening processes during an elasto- 
plastic deformation of materials. 

Consider, for instance a homogeneous isotropic case. Let Qu be a symmetric 
(0,2)-tensor. The evolution in the direction of this tensor, i.e. the evolution for 
which Efj = X(t)Qu, X(t) > may proceed only if the difference (Tfj - Tf}) is 
such that Tr g ((Tfj — T}j)Q IJ ) ^ 0, i.e. if this difference is positive in the direction 
of tensor Qij. 

Leaving further study and comparison of these conditions with the usual yield 
criteria, (36] 03] for future work, we notice only that the conditions (6.18) are 
anisotropic by its nature and might possibly provide useful supplement to the usual 
criteria in essentially anisotropic situations. 

7. Conclusions 

In this work we analyzed the relation between the Bilby-Kroner -Lee multiplica- 
tive decomposition F = F e F p of the total deformation gradient into elastic and 
plastic factors ( [26l H2 OH HB] ) and the theory of uniform materials ( [40 ] 133 ] HT] ) . 
We prove that the Bilby-Kroner-Lee multiplicative decomposition is equivalent to 
the uniform material model with two deformation mappings, i.e. the total <p and the 
inelastic <f>\ deformations together with the uniformity structure. Uniformity enters 
through the (l,l)-tensor field D in the material manifold M or through the mate- 
rial metric g. We introduced the total, the elastic and the inelastic strain tensors 
characterizing different types of the geometrical evolution of the material. After 
discussing the relations between these strain tensors and the deformation gradients 
F e and F p , we chose the form of the internal energy (5.2) and of the dissipative 
potential (5.6) for the materials modeled by the triple ((f), <^i,g). The evolution 
equations were written down for all dynamical variables (4>,4>i,g). We discussed 
different types of stress tensors that naturally enter the scheme of our work. Fi- 
nally we wrote down the dissipative inequalities for the materials of (</>, (j>i , g)-type 
where the terms corresponding to the different types of dissipative processes are 
separated. 

Further research along the lines indicated in this paper seems to be in order. 
First, one should compare our results with those obtained by G.Maugin in a different 
framework [25] 124] . Second, in the continuation of this work we will study the 
evolution equations (5.7-5.12), obtain the energy balance law and the heat equation 
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that follows from it along the lines of [25 . Third, some special cases and examples 
will be considered. 

Acknowledgements. The authors would like to thank professor M.Elzanowski 
who made the manuscript of the book [TT] available to them before the publication 
and for valuable advices during the preparation of the manuscript. We also 

8. Appendix 

In this Appendix we present calculation of the total strain rate E tot that was 
used in Sec. 6. 

From the formula (4.2) for the linearized definition of strain tensors we get 

E tot = gz X C{<h)E* + g^gE m + E m . 

Taking derivative we get 

E tot = g^Ci^E* 1 + g^C(4> 1 )E a + g^gE™ + g^gE™ + E m . (8.1) 

From the definition of linearized E m — \g~ 1 {C{4>\) — g) we get C(4>i) = .9 + 2gE m 
and, therefore, <7(0i) = g + 2gE in + 2gE m . As a result, 

g^Cifa) = g \g + 2gE m + 2gE m ) = 2E m + 4E m E m + 2(1 + 2E m )E m , 

where we have used g$ g = 1 + 2E m . 

In the second term in (8.1) g^Cfa) = g^gg^C^i) = (l+2E m )(l+2E in ), in 
the third one g^ 1 ^ — 2E m . Substituting these expressions into (8.1) and collecting 
coefficients of strain rate tensors we get 

E tot = (1 + 2E m ) (1 + 2E in )E el + (1 + 2E m )E m {\ + 2E el ) +E m (l + 2E el )(l + 2E in ). 

(8.2) 

In a case where strain tensors participating in the second and third terms of the 
last formula commute with the corresponding strain rate tensor, we get 

E tot = (1 + 2E m )(l + 2E m )E el + (1 + 2E m ){\ + 2E el )E in + (1 + 2E el )(l + 2E in )E m 

(8.3) 
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